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Abstract — We consider a single-antenna Gaussian multiple- 
access channel (MAC) with a multiple-antenna amplify-and- 
forward (AF) relay, where, contrary to many previous works, also 
the direct links between transmitters and receiver are taken into 
account. For this channel, we investigate two transmit schemes: 
Sending and relaying all signals jointly or using a time-division 
multiple-access (TDMA) structure, where only one transmitter 
uses the channel at a time. While the optimal relaying matrices 
and time slot durations are found for the latter scheme, we 
provide upper and lower bounds on the achievable sum-rate 
for the former one. These bounds are evaluated by Monte Carlo 
simulations, where it turns out that they are very close to each 
other. Moreover, these bounds are compared to the sum-rates 
achieved by the TDMA scheme. For the asymptotic case of high 
available transmit power at the relay, an analytic expression is 
given, which allows to determine the superior scheme. 

I. Introduction 

In today's wireless communication systems, the demand for 
higher data rates and wide-range coverage is steadily growing. 
To meet these requirements, a high density of base stations 
is necessary, which entails high costs for installation and 
maintenance. Another possibility to increase throughput and 
coverage is the use of relay nodes, which have much lower 
costs. Relay channels were considered in [1] first, and have 
drawn more and more research attention in the last decades. 

Depending on how the signals are processed at the relay, 
different types of relaying schemes are distinguished. The most 
common ones are amplify-and-forward (AF, also called non- 
regenerative relaying) and decode-and-forward (DF, also called 
regenerative relaying). While in AF, the relay simply amplifies 
the received signals subject to a power constraint, a complete 
decoding and re-encoding of the signal is necessary when 
using DF. As this yields higher costs and larger delays, we 
will restrict ourselves to AF relaying schemes in this paper. 

For multiple input multiple output (MIMO) systems with 
additive white Gaussian noise (AWGN), the main challenge is 
to find both the covariance matrix at the transmitter and the 
matrix that maps the relay's input to its outputs, such that the 
data rate is maximized. The problem becomes even harder 
to solve, if a relay system with multiple transmitters, also 
called a multiple-access relay channel (MARC), is considered. 
This holds especially if the direct links between transmitters 
and receiver are also taken into account. A solution for this 
general problem has not been found yet. However, numerous 
previous works have made considerable progress at least for 
some simplified versions of the problem: 



If the direct links are neglected, the optimal structure of 
both the relaying matrix and the transmit covariance matrices 
has been found (21 0. However, this structure still contains 
parameters that are subject to optimization and the optimal 
solution of this problem remains unknown. For the case of 
a single receive antenna, the above problem could be solved 
in where it was also shown that time-division multiple- 
access (TDMA) further increases the achievable sum-rate. In 
0, a single-user system was considered, where the transmit 
covariance matrices were fixed to scaled identity matrices. 
With this restriction, an algorithm was found that optimizes the 
relay matrix. However, for the case of non-zero direct links, 
only upper and lower bounds could be provided. Different 
from all previously mentioned works, a half-duplex relay was 
assumed in (6). This relay was used in single-user systems 
both with and without direct links, for which suboptimal 
transmit strategies based on iterative algorithms were derived. 

In this work, we consider a full-duplex if -user MARC 
with an AF-relay and non-zero direct links, where only the 
relay has multiple antennas. For this system, we first derive 
new upper and lower bounds on the achievable sum-rate 
for the case where all transmitters send their signals jointly. 
Subsequently, we will extend the TDMA-based transmission 
scheme introduced in (4) to the case where the direct links 
are present. An optimal solution for this scheme is achieved 
by an iterative algorithm. Finally, the achievable sum-rates of 
the TDMA and the "joint relaying" scheme will be compared, 
where it can be seen that the superiority of TDMA found in |4 1 
does not always persist for the case of non-zero direct links. 

This paper is structured as follows: In Section [II] we intro- 
duce the channel model and describe the constraints that have 
to be fulfilled while optimizing the sum-rate. Subsequently, 
we derive upper and lower bounds for the joint relaying 
scheme in Section [III] The TDMA scheme is discussed in 
Section IIVI where also an algorithm achieving the optimal 
solution and an asymptotic comparison to joint relaying is 
given. Further comparisons for more general scenarios are 
provided in Section[V]by means of simulation results. Finally, 
Section [VTl concludes the paper. 

II. Channel model 

A. Notation 

We denote all column vectors in bold lower case and 
matrices in bold upper case letters. The trace, determinant, 
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Fig. 1. A"-user multiple-access relay channel (MARC) 

Hermitian, and transpose of a matrix A are identified by tr(A), 
|A|, A H , and A T , respectively. We use ||x|| to denote the 
Euclidean norm of a vector x and I to describe the identity 
matrix. Furthermore, A max (A) and v max (A) indicate the largest 
eigenvalue of a matrix A and its corresponding eigenvector. 

B. Channel Model 

The MARC that we consider is depicted in Figure Q] and 
consists of K transmitting nodes. Each user k £ {1, . . . , K} 
transmits the signal x^ k ' £ C, which reaches both the relay 
and the receiver. The channel matrix for the transmission to 
the relay is given by the vector h£ £ C Mr , while the channel 
to the receiver is described by the scalar £ C. Thus, the 



received signal y,. of the relay can be written as 



k=l 



where z r ~ CAT (0,1) is the additive white Gaussian noise at 
the relay and M r denotes the number of antennas at the relay. 
The relay amplifies the signals by the matrix F and transmits 
the signal x r = Fy r over the channel h H £ C lxMr to the 
receiver. 

It is assumed that the transmission from the relay to the 
receiver takes place in a different frequency band, i.e., the 
signals transmitted by the relay are orthogonal to the signals 
transmitted by the users. To overcome the problem that the 
signals from the relay arrive with the delay of one symbol, we 
assume that the direct signal can be buffered. Hence, the out- 
of-band reception can be modeled by a virtual second antenna 
at the receiver, which receives only the signal from the relay. 
This results in the received signal 



y = 



K 

h H F £ hrV fe ) + h H Fz r + zi 

k=l 



K 

E 

k=l 



where Zi ~ CJ\f(0,Ng) (i — 1,2) denote the Gaussian noise 
terms at the receiveiQ As the first component of y contains 
noise both from the relay and from the receiver, it can be 

'Throughout this paper, we assume TVq = 1 



simplified by normalization without changing the systems 
properties. Thus, we will use 

K 

(1) 



/c=i 



as equivalent output at the receiver, where z ~ CJV(0, NqT), 

h[ k) = [r-^h^Fh^, hfY, and r = 1 + h^FF^h. 

Both the transmitters and the relay are subject to average 
power constraints, which are given by 

,, , l|2 x 



< pW Vke {I,..., K} 

E(tt (x r xf )) =tr (f U + J2 h< fc) P (fe) h«^ F*j < P r . 

(2) 

Finally, we assume that perfect channel state information is 
available at all nodes. 

III. Joint Relaying Scheme 

The transmit strategy, where all transmitters send their 
signals at the same time will be referred to as "joint relaying" 
in the remainder of this paper. In this case the MARC from 
([TJ can be interpreted as a pure single input multiple output 
(SIMO) multiple-access channel (MAC). For this SIMO MAC, 
the achievable sum-rate can be optimized by influencing the 
channel gain through the choice of the relaying matrix F. As 
in the MAC, the sum-rate can be calculated as 



i? s = log 2 



A' 



H 



k=l 



(k) 

Evaluating this determinant by using the definition of h eS 
and standard transformations of linear algebra, i?s can be 
reformulated as 



i?s = log 2 (1 + s + r-'h^F ((1 + s) R - T) F^h) 
= log 2 (l + s + r^h^F (R + W) F H h) , 



(3) 
(4) 



where s 
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P {k \ R = Ef =1 hl fc) P^h [ r k)H , and 



\k=l , 
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\k=l J 



H 



(*0) pO')p(fe). 



j,k=l 

As it can be seen from the above equations the choice of F 
only influences the last term inside the logarithm in (0 and 
©, while the other terms are constant. However, compared 
to the optimization problem with absent direct links 0), we 
have the additional term W, which occurs in the sum-rate but 
not in the power constraint (0. Hence, the optimal relaying 
matrix is not the same as in [4|. As the optimal solution seems 
to be hard to find, we will derive upper and lower bounds in 
the following two subsections. 



A. Upper Bounds on i?s 

A first upper bound can be obtained from the fact that T is a 
positive semidefinite matrix. Thus, i?s can be upper bounded 
by 



Rs < log 2 (l + s + r~ 1 h H F ((1 + s) R) F H h) 



(5) 



Besides the additive term s and the constant factor (1 + s) the 
maximization problem is now similar to the one in [4], such 
that this upper bound is optimized by choosing (cf. J2] 0) 



P 



(6) 



yi + A max (R) ||h 

Using this relaying matrix in the right side of (0, we obtain 

P W x.JlJ A max (R)|lh|| 2 P r 
i?S,up,l = log 2 (1 + S) 1 + 2 — 

y y 1 + ||h|| P r + A max (R) 

A second upper bound can be obtained by ignoring the relay 
power constraint, i.e., by letting P r — > oo, which delivers 

h H F (R + W) F H h" 



i?S < l0g 2 [1 + S+ h H FF H h 

< log 2 (l + s + A max (R + W)) = (7) 

where the second inequality follows from the Rayleigh quo- 
tient. By design, -Re, U p,2 becomes tight at high values of P r , 
while i?s,up,i is tighter if P r is small. 

B. Lower Bounds on i?s 

In order to find rates that are actually achievable, it is 
possible to choose F as in the derivation of the upper bounds, 
although these choices will not be optimal in general. One 
possibility is to choose F as in ([§). Another approach, which 
is derived from the second upper bound i?s jUp ,2, is to set 



F = 7-^-v^ ax (R + W), 



where 7 £ K is chosen such that the relay power constraint 
is fulfilled with equality. Throughout all numerical simulations 
that have been made, the second approach ([HJ turned out 
to deliver better results. For this reason, it will be the only 
considered lower bound in the remainder of this work. If F is 
chosen as in ((HJ, the achievable rate can be written as 



ife = log 2 1 + s + 



|h|r A max (R + W) 7 2 
l + l|h||V 



(9) 



IV. TDMA-Based Relaying 



In this section, we will introduce a relaying scheme based 
on TDMA as in H. This scheme includes a division of the 
transmission in K time slots, where user k occupies the k- 
th time slot exclusively. Also the relay incorporates this slot 
structure, i.e., the relaying matrix F^ in time slot k can be 
adapted to the channel of user k only. Thus, the TDMA slot 
structure decomposes the channel in K independent single- 
user relay channels. Therefore, we will first derive the optimal 
structure of the transmit covariance- and relaying matrix for 



the single-user relay channel in Subsection IIV-AI In the 
following Subsection IIV-BI we will transfer this scheme to 
the MARC with TDMA and derive an algorithm that finds the 
optimal duration of the time slots, such that the sum-rate is 
maximized. Finally, in Subsection IIV-CI we will compare this 
sum-rate to those derived in section [Til] for the case P r — > 00. 

A. Single-User Relaying 

In order to describe a single-user relay channel with a 
consistent notation, we assume the same channel model as 
introduced in Subsection IH-BI with only K = 1 transmitting 
user. Thus, also (0]) is valid and can be used to calculate the 
(sum) rate R^ 1 ' of the only user. In contrast to the previous 
section the optimization of F^ is strongly simplified as we 



have R = h { r 1] P^h 



H 



p( x \ and especially 



W = 0. Hence, the rate can be written as 



h^FWRFW^h 
1 + h«F( 1 )F( 1 ) ff h 



= log 2 1 + s + 



(10) 



The optimization of (TTOb over F^ is basically the same as in 
(0. Thus, in analogy to (0 the optimal F' 1 ' is given by 



p(i) 
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which leads to a rate of 
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RW = log. 
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p(i). 
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p(i) 



(8) B. TDMA-based Transmission Scheme 



The If -user MARC is decomposed in K single-user relay 
channels by using a TDMA scheme, such that user k transmits 



only in a time slot of duration r' fe ^ > with Yl,k=i = ^■ 
In each time slot, the optimal choice of the relay matrix F^ fc ^ 
can be obtained as in subsection IIV-AI The only difference is 
the transmit power constraint: As user k only transmits in 
fraction of the time, it can use a transmit power of p( fc )/r^ fe ^ 
and still fulfills the average transmit power constraint (0. 
Thus, the rate of user k is given by 



R (k) _ T (k) b l + 



p(k) 



|h|| 2 p( fe ) 



+ 



(k) 



P 



Ihll 2 P.t« +P( fc ) h. 



Ak) 



(11) 



Ak) 



and the sum-rate can be calculated as Psjdma = 
X}fc=i-^ fe ^ This sum-rate can be optimized by the choice 



of , . . . , t^ k \ i.e., we are facing the optimization problem 

K 



.(fe) 



S.t. %)=l-^T< fe ' =0, 



(12) 



fc=l 



where r = [r^ 

Vj 7^ A; and 



(K) 



] . It is easy to see that 



on 



(fc) 



drU) 



= 



OtW 



< 0, which makes the problem con- 
vex. Thus, the famous Karush-Kuhn-Tucker (KKT) conditions 
provide necessary and sufficient conditions for optimality. 
For the above problem, the KKT conditions of a solution 
t* to be optimal can be formulated as H(t*) = and 
Vi? E ,TDMA(T*) + v*Wh(r*) = 0, where v* E K can be 
chosen arbitrarily. As the derivatives of h are directly obtained 
as 



Oh 
0r< k ) 



-1, the second KKT condition can be rewritten as 

' Vfc, 



<9t« 

i.e., the derivatives of the individual rates have to be the 
same for each user. Due to their lengthiness, those derivatives 
are not stated here. However, they are straightforward to 

8flW ' 



calculate and it is easy to see that 

3fi< fc > 



T C*)=0 



oo and 



= 0. Unfortunately, a closed form solution how 

— ► oo 

to optimally choose r' 1 ), . . . , t^ k ' as in J4| seems intractable. 
Therefore, we describe in Algorithm Q] how the optimal 
solution can be found iteratively. 



Algorithm 1 Iterative optimization of . . . , r 



1: Set r«* = i Vfc = l, 
while true do 

on (k) 

i = arg mm fc 
j = argmaxfe 

' t W=tW* 



4: 

5: 
6: 



<fc) =r (fc)* 



if 



OtU) 



(k)— T {k)» 
on^ 



> e then 

t W)=tW)* r< i )=r( i )* 

Find < <5 < min{rW, 1 - r®}, such that 



7: 
8: 
9 

10: 
11 
12 



9r«) r ( 3)=T ( 3 ). +(5 

T w = T w*_ s 

else 

break 
end if 
end while 



The main idea of the algorithm is to iteratively equalize 
the derivatives of by changing the lengths of the time 
slots tW, . . . ,t( k '. Therefore, the users i and j with the 
smallest and largest derivative are selected. Their derivatives 
are equalized by numerically finding a value 6, which is added 
to and subtracted from r^. Due to the properties of the 



functions discussed above this value can always be found in 
the interval (0, min-fr^. 1 — r' 3 '}), and the other derivatives 
remain unchanged. This procedure is repeated iteratively until 
the difference of the largest and smallest derivative is at most 
e, which can be selected very small to approximate the optimal 
solution as good as desired. 

C. Comparison with Joint Relaying 

As a closed form of the optimum sum-rate of the TDMA 
scheme has not been found, and for joint relaying only upper 
and lower bounds are available, comparing these two schemes 
in an analytical way is not as straightforward as in H. 
Therefore, the two schemes will mainly be compared by the 
mean of simulation results in Section [V] However, an analytic 
comparison is possible for the asymptotic case of P r — > oo, 
which is described in the following theorem. 

Theorem 1: In the considered X-user MARC with direct 
links and P r — > oo, the joint relaying scheme achieves higher 
sum-rates than the TDMA-based transmission scheme with 
optimal time slot durations T^ k \ if and only if 



A max (R + W) > 



K 

E 

fc=i 



p(fc)_ 



(13) 



Proof: For the joint relaying scheme with P r oo, we 
have 7 — > oo in i.e., the sum-rate Rs converges to -Rs, U p,2 
in (O. Considering the TDMA scheme for P r — > oo, the 
individual user rates R( k > from ( fTTT i tend to 

p(k) 



Rt ] = r« lo. 



'82 



1 



(k) 



If this term is used instead of R^ in the optimization problem 
( fT2l . it is straightforward to show that choosing r^ k ' as 



p(k) 



K 

g Pti) 
3 = 1 



Aj) 



is the global optimal solution that leads to the sum-rate 

K 



Re 



TDMA, 00 



= log 2 



1 



fc=i 



(fe) 



(14) 

Comparing ( fl4l > and ((7]), the theorem follows. ■ 
Considering especially the matrix W in $13[ , it turns out that, 

(k) 

if the direct links h\ are zero, we have W = 0. Thus, the 
left hand side in (TPTl i is smaller and TDMA is always better in 
this case (cf. [4|). However, if the direct links are increased, 
also A max (R + W) increases rapidly, which compensates the 
disadvantage of joint relaying. 

V. Simulation Results 

In order to see how the strength of the direct links affects the 
superiority of either TDMA or joint relaying, we will evaluate 
the achievable sum-rates by Monte Carlo simulations in this 
section. Moreover, the average gap between upper bounds and 
achievable rates shall be investigated. 
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Fig. 2. Comparison of TDMA-based and joint relaying scheme for K = 10, 
M r = 4, P max /7Vo = 10 dB 
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Fig. 3. Probabilities of joint relaying being better than TDMA at P r — > oo 
for K = 10, M r = 4 

For this purpose, we assume channels with independent 
Rayleigh-fading, i.e., all entries of the channel matrices 
are ~ CAf(0, 1) and independent from each other. As an 
exception, the channel gains for the direct links h^' are 
multiplied with a parameter aet, i.e., they are ~ CAf(0, a 2 ) 
in order to vary the strength of the direct links. All results 
presented in this section are obtained by averaging over at 
least 1000 channel realizations and the transmit powers P^ 
of the users are assumed to be uniformly distributed between 
and P max . 

For a system with K = 10 users, the sum-rates of both 
TDMA and joint relaying are plotted in Figure [2] We assumed 
that the relay has M r = 4 antennas and that the maximum 
transmit power to noise ratio is P max /No = 10 dB, while 
both the strength of the direct channels a and the power at 
the relay P r are varied. It can be directly observed that for 
all parameters the minimum of the two upper bounds and the 
achievable rates of the joint relaying scheme are very close 
to each other. Concerning the comparison with TDMA, it can 



be seen that for a = 0.1, the superiority of TDMA discussed 
in [4 1 persists only for large values of P r . If the strength of 
the direct links is further increased, the sum-rates achieved by 
TDMA grow only slowly, especially for large values of P r . 
On the other hand, the sum-rates of the joint relaying scheme 
grow faster, such that they are clearly higher than those of 
TDMA for a = 0.3 and especially for a = 1. 

A further comparison of TDMA and joint relaying for 
different maximum transmit powers is given in Figure[3] In this 
figure, we plotted the probability that joint relaying achieves 
higher sum-rates than TDMA at P r —> oo, i.e., the probability 
that ( TT3l > holds, for different values of a and P max . As already 
observed in Figure [2] these probabilities grow with the strength 
of the direct links a. Moreover, it can be noticed that the 
probabilities also grow with increasing transmit powers. Note 
that, also if dT3l > does not hold, joint relaying can still achieve 
higher rates for lower P r (cf. plots for a = 0.1 in Figure |2). 

VI. Conclusion 

We have considered a K-usei MARC with direct links, 
where only the relay has multiple antennas. For this channel, 
we found upper and lower bounds on the achievable sum-rate 
when all stations transmit their signals jointly. It was shown by 
Monte Carlo simulations that these bounds can be very close 
to each other. If a TDMA protocol is used, we were able 
to determine the optimal relaying matrices and obtained 
the optimal time slot durations by an iterative algorithm. For 
the asymptotic case of P r — > oo, an analytic expression was 
derived, which allows to determine whether joint relaying or 
TDMA achieve a higher sum-rate. For general values of P r 
the performance of joint relaying and TDMA were also found 
by Monte Carlo simulations. It turned out that the superiority 
of TDMA found in |4| for absent direct links gets lost quickly 
if the strength of the direct links is increased. This trend is 
accelerated if the available transmit powers are increased. 

References 

[1] T. Cover and A. E. Gamal, "Capacity Theorems for the 
Relay Channel," Information Theory, IEEE Transactions on, 
vol. 25, no. 5, pp. 572-584, Sep. 1979. [Online]. Available: 
http://ieeexplore.ieee. org/xpls/abs_all.jsp?arnumber= 1056084 

[2] Z. Fang, Y. Hua, and J. C. Koshy, "Joint Source and Relay Optimization 
for a Non-Regenerative MIMO Relay," in Sensor Array and Multichannel 
Processing, 2006. Fourth IEEE Workshop on. IEEE, Jul. 2006, pp. 239- 
243. [Online]. Available: http://dx.doi.org/10.1109/SAM.2006.1706129 

[3] Y. Yu and Y. Hua, "Power Allocation for a MIMO Relay System 
With Multiple-Antenna Users," IEEE Transactions on Signal Processing, 
vol. 58, no. 5, pp. 2823-2835, May 2010. [Online]. Available: 
|http://dx.doi.org/ 10.1109/TSP.2010.2042476l 

[4] F. Rnabe, O. Mohamed, and C. Huppert, "Superiority of TDMA 
in a Class of Gaussian Multiple-Access Channels with a MIMO- 
AF-Re lay," arXiv:1202.1734 v2, Mar. 2012. [Online]. Available: 
http://arxiv.org/abs/1202.1734 

[5] X. Tang and Y. Hua, "Optimal Design of Non-Regenerative MIMO 
Wireless Relays," Wireless Communications, IEEE Transactions on, 
vol. 6, no. 4, pp. 1398-1407, Apr. 2 007. [Online]. Available: 
http://dx.doi.org/10.1109/TWC.2007.348336 

[6] N. Varanese, D~ Simeone, Y. Bar-Ness, and U. Spagnolini, 
"Achievable Rates of Multi-Hop and Cooperative MIMO Amplify- 
and-Forward Relay Systems with Full CSI," in Signal Processing 
Advances in Wireless Communications, 2006. SPAWC '06. IEEE 
7th Workshop on. IEEE, Jul. 2006, pp. 1 -5. [Online]. Available: 
http://dx.doi.org/10.1109/SPAWC.2006.346423 



